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Resonances and bifurcations in axisymmetric scale-free 
potentials 
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ABSTRACT 

We investigate an analytical treatment of bifurcations of families of resonant 'thin' 
tubes in axisymmetric galactic potentials. We verify that the most relevant bifurcations 
are due to the (1:1) resonance producing the 'inclined' orbits through two different 
mechanisms: from the disk orbit and from the 'thin' tube associated to the vertical 
oscillation. The closest resonances occurring after these are the (4:3) resonance in 
the oblate case and the (2:1) resonance in the prolate case. The (1:1) resonances are 
treated in a straightforward way using a 2nd-order truncated normal form. The higher- 
order resonances are instead cumbersome to investigate, because the normal form has 
to be truncated to a high degree and the number of terms grows very rapidly. We 
therefore adopt a further simplification giving analytic formulas for the values of the 
parameters at which bifurcations ensue and compare them with selected numerical 
results. Thanks to the asymptotic nature of the series involved, the predictions are 
reliable well beyond the convergence radius of the original series. 

Key words: galaxies: kinematics and dynamics - methods: analytical. 



1 INTRODUCTION 



, Axisymmetric galactic po tentials admit only tube orbits 
■ around the symmetry axis (iBinney fc Tremainellioosi ') . How- 
ever, departing from the two linear approximations around 
, the circular orbit, epicyclic motion in the symmetry plane 
' and the oscillation perpendicular to this plane, nonlinear 
dynamics are characterized by a hierarchy of resonances. 
To these are associated the bifurcations of different kinds 
of 't hin' tubes which p arent the corresponding 'thick' fami- 
lies l|Contopoulosll2004l ). The occurrence of this phenomenon 
may heavily affect the kinematics of stars with, for exam- 
ple, sudden removal from the disk (levitation) and/or veloc- 
ity redistribution that can be viewed as a heating mecha- 
nism of the disk l|Sridhar fc Touma|[l996h . Resonant fami- 
lies may also play a relevant role in self-consistent models of 
spheroidal halo or disk/halo systems. 

Scale-free potentials arc often adopte d as realistic galac- 
tic models tRich stone 198 2; Eva ns 1994: Touma fc Tremain^ 
Il997l : ISridhar fc Toumal Il999l ) The orbit structure is the 
same at each energy level, resulting in a much simpler 
description of dynamics. Their axisymmetric version is 
parametrized by only two parameters: a, the exponent of 
the power law (which can be taken zero in the limit case 
of the logarithmic potential) and q, the ellipticity of the 
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isopotentials. Adding to these L, the conserved angular mo- 
mentum around the symmetry axis, completely specifies the 
orbit structure. The ensuing 2-dof dynamical system is in 
general non-integrable. An analytic approach can be at- 
tempted by expanding the effective potential around the 
critical point corresponding to the circular orbit and push- 
ing the expansion beyond quadratic terms. The Hamilto- 
nian series is afterwards replaced by a 'normal form' which 
plays the role of an integrable appro ximat i on of the original 



non-integrable system (|Contopo ulo3 19631 : iGustavsMilllQGd : 
iGerhard fc Sahalll99ll : mmguasil200ll 'l. The dynamics of the 



normal form are amenable to a totally analytic treatment 
that allows us, among other things, to find the bifurcation 
thresholds of periodic orbits in terms of the param eters of 
the effective potential (jBelmonte et al.ll2006l . [20o3) and ap- 
proxima te expressions of the solutions of the equations of 
motion l|Pucacco et al.ll2008bl '). 

In the present work, we investigate the properties of the 
simplest non-trivial series expansion of axisymmetric scale- 
free potentials beyond the epicyclic approximation. We ver- 
ify that the most relevant bifurcations are due to the (1:1) 
resonance producing the 'inclined' orbits through two dif- 
ferent mechanisms (Hunter et al,.1998. ): from the disk orbit 
(the 'horizontal' normal mode) and from the 'thin' tube as- 
sociated to the z oscillation (the 'vertical' normal mode) 
The closest resonances occurring after these are the (4:3) 
resonance in the oblate case and the (2:1) resonance in the 
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prolate case. Other high-order resonances can be present at 
very low values of the angular momentum. 

The (1:1) resonances are treated in a straightforward 
way using a 2nd-order truncated normal form. The higher- 
order resonances are instead cumbersome to investigate, be- 
cause the normal form has to be truncated to a high degree 
and the number of terms grows very rapidly. We therefore 
ad opt a further simplification suggested by the approach 
of ISanders et all |2007, 1. We give therefore analytic formu- 
las for the values of a, q, L at which bifurcations ensue and 
compare them with selected numerical resu lts. Thanks to 
the as ymptotic nature of the series involved l|Pucacco et al.l 
l2008al 'l. the predictions are reliable well beyond the conver- 
gence radius of the original series. 

The plan of the paper is as follows: in section 2 we 
describe the family of systems under investigation and how 
they are prepared to undergo a perturbative analysis; in sec- 
tion 3 we recall the method based on Hamiltonian normal 
forms; in section 4 we apply the theory to compute bifur- 
cation thresholds of the main resonant families; in section 5 
we conclude with a discussion of the results. 



2 THE POTENTIALS AND THEIR SERIES 
EXPANSIONS 

We investigate the dynamics of the family of potentials 

a/2 



(1) 



l\og[R^+-^y a = 0. 

The ellipticity of the equipotentials is determined by the 
parameter q: we have an 'oblate' figure when q < 1 and a 
'prolate' figure when q> 1. The slope of the power law will 
be restricted to the range 

- 1 < a ^ 2. (2) 
The Hamiltonian of the system in cilindrical coordinates is 



(3) 



that, exploiting the conservation of the axial angular mo- 
mentum 

P<t> = L, (4) 

is effectively the Hamiltonian of a two degrees of freedom 
system in the family of potentials 



L 

Vc{R,z-L,q) = ^ -f $Q(ii,z;g). 



(5) 



These potentials have a unique absolute minimum in 

R = Rc{a)=L^, 2 = 0. (6) 

This is a stable equilibrium on the meridional plane (j> — 
const corresponding to a circular orbit of radius Rc{a) of 
the full three-dimensional problem. Since the dynamics are 
scale-free we may fix energy once and for all 



E — Ea = 



(i+i)e-— , «/0, 
0, a = 0. 



(7) 



This implies that the radius of the circular orbit at this 
energy is 



Rc{a) = e"~, -1 < a < 2 

and we can investigate the dynamics at 

E = Ea, Vae (-1,2], 

by varying L in the range 

1 



0< L^Ln 



Rc 



(8) 
(9) 

(10) 



without any loss of generality l|Hunter et al.|[l998l ). 

In order to implement the perturbation method we ex- 
pand the effective potential around the minimum A 
common attitude in normal form theory holds up the relia- 
bility of its predictions well beyond the convergence radius 
of the expansion. Although a rigorous proof of this state- 
ment is still lacking, it is confirmed by several resul ts ob- 
tained in analytical and celestial mechanics (Scufiair^ ll995l : 
lEfthvmiopoulos et all l2004l : IPucacco et~all 2008a). One of 
the aims of the present work is to test the reliability of these 
asymptotic estimates. 

We introduce rescaled coordinates according to 

. R — Rc .2 . , 

y^-FT (11) 



Rc 



with origin in the equilibrium point (|6]). The potential ([Sjl is 
expanded as a truncated series (in the coordinates x, y) of 
the form 



V^^\x,y;L,q) = Cjj.k-j) (L, a, q)x^y''~ 



(12) 



fe=0 j=0 



where the truncation order TV is determined by the reso- 
nance under study and is discussed below. From © and the 
rescaling (|11[) . the constant term of the expansion is 



C(o,o)iL,a) 



(1 + ^)^*^. «/o, 
i +L, a = 



and the other coefficients have the form 

2a 

C(j,k-i){L,a,q) = L2+a C(j_fc_j) (a, g). 



(13) 



(14) 



In order to simplify formulas, we introduce the new param- 
eter 



= - 



2a 

2 + a' 



-l</3^ 2, 



(15) 



with the same range of a in view of ((5)) . 

The orbit structure of the original family of potentials 
((5]) at the energy level fixed by ((7)l will be approximated by 
the orbit structure of the rescaled Hamiltonian 



H =^{pl +pI) + Vi'^\x,y-q), 
where 

V^''\x,y-q)^L^V^''\x,y-L,q). 



(16) 



(17) 



The dynamics given by Hamiltonian (|16p take place in the 
rescaled time 



at the new 'energy' 



(18) 



i5 = {Ea - C(o,o) {L, a)) = - (l - (L/L,„,,)^) . (19) 
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According to (|10p the singular value L = is excluded from 
the analysis implying that the fictitious energy (|19|l is always 
finite and that the expansion around the equilibrium point 
((6| make sense. 

The nonvanishing coefficients of the expansion of V^i^' 
up to order A'^ = 4 are the following: 



C(2,0) = 

C(0,2) = 

C(3,0) = 

C(l,2) = 

C(4,0) = 

C(2,2) = 

C(0,4) = 



2+a 
2 ' 
1 

2^^' 

10 + 3a-a^ 
6 ' 
2- a 

54+lla-6a^+ct^ 
24 

6 — 5a — a*^ 

— ' 



49-^ 



(20) 
(21) 
(22) 
(23) 

(24) 
(25) 
(26) 



The first two of them provide the frequencies of the epicyclic 
motions. Recalling the time rescaling in (|18|l , the radial and 
vertical harmonic frequencies are respectively 



V2 + ^ 



and 



3 NORMAL FORMS 



(27) 



(28) 



Hamiltonian (I16|l is in the form of a power series and is there- 
fore naturally apt to be treated in a perturbative way as a 
non-linear oscillato r system. We construct a 'no rmal form' 
(jSoccaletti fc Pucac co IQOd lSanders et al.ll2007l \ namely a 
new Hamiltonian series which is an integrable approxima- 
tion of the original one, suitably devised to catch its most 
relevant orbital features. 

The normal form is 'non-resonant' when the two 
harmonic frequencies (|27|l and (|28|l are generically non- 
commensurable: in this case we get explicit formulas for 
actions and frequencies of the box orbits parented by the 
radial (disk) and by the vertical (thin tube) orbits. A 'res- 
onant' normal form is instead assembled by assuming from 
the start an integer value for the ratio of the harmonic fre- 
quencies and including in the new Hamiltonian terms de- 
pending on the corresponding resonant combination of the 
angles. This possibility might be considered an exception: it 
is instead the rule because, even if the unperturbed system 
is non-resonant with a certain real value 



p — k/u 



(29) 



of the frequency ratio, the non-linear coupling between the 
degrees of freedom induced by the perturbation, determines 
a 'passage through resonance' with a commensurability ra- 
tio, say mi/m2, corresponding to the local ratio of oscilla- 
tions in the two degrees of freedom. This in turn is respon- 
sible of the birth of new orbit families bifurcating from the 
normal modes or from lower-order resonances. 

To evaluate the most relevant resonances in our case, 
let us come back to the epicyclic frequencies (|27p and 
Their ratio is 



We then approximate the frequencies with a rational number 
plus a small 'detuning' llContopoulos fc Moutsoula3 Il966l : 
Ide Zeeuw fc Merritt|[l98l ) 

P=^+5 (31) 
m2 

and proceed like as the unperturbed harmonic part would 
be in exact mi : m2 resonance putting the remaining part 
inside the 'perturbation'. We speak of a detuned (mi:m2) 
resonance, with 



Nmin = mi -f 7712 



(32) 



the order of the resonance. In general, a normal form trun- 
cated at A'^min includes the first resonant term. With a in 
the range ([2]) and an interval 



0.5 < g < 1.5 



(33) 



(34) 



of reasonable values of g, we see that the most relevant res- 
onance values with low commensurability are 

1 2 3 1 4 3 2 
2' 3' 4' T' 3' 2' T' 
In our analysis, we will mainly concentrate on the 'central' 
value 1 : 1 and on the higher-order cases 2 : 1 and 4 : 3. In 
fact, periodic orbits with these fequency ratios and associ- 
ated quasi-periodic families are usually the most prominent 
in nu merical investigations (|Hunter et al ]| 19981 : IContopoulod 
l2004f ). The procedure can be easily extended to other cases. 

Normal forms for the Hamiltonian system corre- 
sponding to (1161) are const ructed with standard methods 
(Boccaletti fc Pucaccd Il999[ ) and were used to determine 
the main features of the orbit structure of the logarit hmic 
potential (|Belmonte et al.ll2007l : IPucacco et al.l 12008^ 1. We 
briefiy resume the procedure in order to fix notations. After 
a scaling transformation 



X — > \/2 + ax, px — > Px/^fi + a, (35) 

y — Py — * ^/qpy, (36) 

the original Hamiltonian ([16} undergoes a canonical trans- 
formation to new variables Px,Py,X and Y, such that 



K(Px,Py,X, F) =^A-„, 

n = 

generated by a function of the form 
G ~ Gi -j- G2 + ... 

with the prescription {K in 'normal form'' 
{Ko,K} = 0. 



(37) 



(38) 



(39) 



In these and subsequent formulas, we adopt the convention 
of labeling the first term in the expansion with the index 
zero: in general, the 'zero order' terms are quadratic homo- 
geneous polynomials and terms of order n are polynomials 
of degree n -I- 2. The zero order (unperturbed) Hamiltonian, 







Ho = -{uJi{PI+X^)+uj2{Py+Y'')) 



with 'unperturbed' frequencies 

1+/3 



hi L 



(jJ2 



V L 



1/g, 



(40) 



(41) 



p = g\/2 + a. 



(30) 



plays, by means of the fundamental equation (|39p , the dou- 
ble role of determining the specific form of the transforma- 
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tion and assuming the status of the second integral of mo- 
tion. The lowest order term of the generating function, Gi, 
is a cubic polynomial. 

Using 'action-angle'-/«fce variables J, 9 defined through 
the transformation 



X 
Px 



(42) 
(43) 



the typical structure of the resonant normal fo rm (truncated 
when the first resonan t term appears) is (jSanders et al.l 
120071 : IContopoulosll2004l ) 



K 



miJi + 7712 J2 + 



mi+m2 

E 



P''='(Jl,J2) + 



amimaJr^JT^ cos[2(m26'i -771162)], (44) 

where "P'*' are homogeneous polynomials of degree k whose 
coefHcients may depend on 5 and the constant amim2{<l, ct) 
is the only marker of the resonance. In these variables, the 
second integral is 

£ = miJi + 7722 J2 (45) 
and the angles appear only in the resonant combination 

tp — 171261 — 171102 ■ (46) 

For a given resonance, these two statements remain true for 
arbitrary A'' > A'min. Introducing the variable conjugate to 



TZ = T772Jl — mi J2, 



(47) 



the new Hamiltonian can be expressed in the reduced 
form K{TZ,ip;£ ,q,a), that is a family of 1-dof systems 
parametrized by £ (and q,a). 

In the applications below we are interested in the global 
structure of phase-space, but the explicit solution of the 
equations of motion is also of great relevance. For a non- 
resonant normal form, the problem is easily solved: the co- 
efficient amim2 vanishes and K no longer has a term con- 
taining angles. Therefore, the J are 'true' conserved actions 
and the solutions are 

X(r) = \/2jrcosKir, Y (t) ^ V2J^ cos{k2T + 9o) , (48) 
where 



(49) 



is the frequency vector and Oq is a suitable phase shift. 

In the resonant case instead, it is not possible to write 
the solutions in closed form. It is true that the dynamics 
described by the 1-dof Hamiltonian K{TZ, ip) are always in- 
tegrable, but, in general, the solutions cannot be written in 
terms of elementary functions. However, solutions can still 
be written down in the case of the 771am periodic orbits, for 
which J, are true action-angle variables. There are two 
types of periodic orbits that can be easily identified: 

(i) The normal modes for which one of the J vanishes. 

(ii) The periodic orbits in general position characterized 
by a /ixed relation between the two angles, 7/7261—777162 = do- 

In both cases, it is straightforward to check that the solu- 
tions retain a form analogous to Ea. (|48|) with known ex- 
pressions of the actions and frequencies in terms of £ and 



the parameters q and a such that «:i/k2 = mi/m2. By us- 
ing the generating function Eq. (|38|) . the solutions in terms 
of standard 'physical' coordinates can be recovered (apart 
from possible scaling facto rs) inverting the canon ical trans- 
formation. As discussed in lPucacco et al.l (|2008lj) the trans- 
formation back to the physical coordinates is expressed as a 
series of the form 



x(t) = Xl + X2 + X3 + ... 

and is given explicitly by 

xi = X, 

X2 = {Gl,X}, 

X3 = {G2,X} + i{Gi,{Gi,X}} 



(50) 

(51) 
(52) 
(53) 



and so forth. From a knowledge of the normalized solutions 
Ea. (|48|l . we can therefore construct power series approxi- 
mate solutions of the equations of motion of the original 
system 

^^-VM''\x,y-q). (54) 

As a rule, normal modes exist on each 'energy' surface 
K = E. Periodic orbits in general position exist instead only 
beyond a certain threshold and we speak of a bifurcation 
ensuing from a detuned resonance. The bifurcation is usually 
described by a series expansion of the form 



E 



E 



Ok • 



(55) 



where the are coefficients depending on the order A'min 
and the parameters q, a. The order of truncation of the se- 
ries is itself related to that of truncation of the normal form 
(|Pucacco et al.|[2008al ). Eq.((55} implies that at exact reso- 
nance (vanishing detuning) the bifurcation is intrinsic in the 
system and that, going away from the initial exact ratio of 
unperturbed frequencies, gradually increases the threshold 
value for the bifurcation. We will see that already a linear 
relation given by the first order truncation provides a reli- 
able estimate and will examine some example of a second 
order truncation. 



4 APPLICATIONS 

We apply the theory resumed in the previous section to a set 
of typical bifurcation problems related to the potentials of 
section 2. In the following we adopt a 'taxonomy' of the main 
resonant families that is coherent with ou r previous work 
and other standard references in the field (|Belmonte et al.l 
I2OO7I : iMiralda-Escude fc Schwarzschildl llOSgl ). In case, we 
mention alternative denominations adopted by other au- 
thors. 



4.1 Non resonant box orbits 

We recall that we are studying the dynamics given by Hamil- 
tonian (|16[l at 'energy' (|19[l : a small value of it corresponds 
to a value of the angular momentum close to its maximum 
amount. This is the quasi harmonic regime in which the 
dynamics in the x — y plane are characterized by box orbits 
oscillating around the two normal modes, Ji = and J2 — 0. 
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The first one is tlie y-axis periodic orbit and, coming back 
to the unsealed variables, it is the 'vertical' z oscillation: in 
3 dimensions it gives the thin tube orbits. The second one 
is the a:-axis periodic orbit that in the unsealed variables is 
the horizontal or 'equatorial' R oscillation and gives the disk 
orbits in 3 dimensions. 

A preliminary step is that of exploiting the non-resonant 
normal form. It can be readily used to get the post-epicyclic 
solution; more interestingly, we will exploit it later to at- 
tempt an approximate description of the higher-order reso- 
nances. The non-resonant normal form of the Hamiltonian 
up to the second order, takes the following form 



K — LO\J\ + UJ2J2 + aJi + hjl -I- cJi J2, 
with 

" " ^ (2C(4,0)'-'? - 5C(3,o)) , 

2cn.2)^i - (3c(i,2)C(3.o) - C(2,2)wf)(tJ? - 4aj|) 

C = ' ■ 



(56) 



(57) 



(58) 



, (59) 



where it is appreciable the appearance of 'small' denomina- 
tors related to the 2:1 resonance, which is the first to appear 
in view of the lowest order coupling term xy^. We also give 
the first term of the generating function psp 



Gi 



C(3,o)Px(2Pl+3X^) 

Wi^^W2(wf — 4a;|) 
2c(i.2)XyfV 



(60) 

(61) 
(62) 



Using expressions (|22H26|I and (|4ip the coefficients in the 
normal form can be written explicitly in terms of the pa- 
rameters of the potential 



-22 + 13a - 

24 ' 
(2-Q)(10-FQ-3g^(2 + «)) 
4g2(2-Fa)(-4 + g2(2 + Q)) ' 
(2-Q)(6-Q-g^(2 + a)) 
2g^/2T^(-4-|-g2(2 + a))• 



(63) 
(64) 

(65) 



The second order post-epicyclic solution is then given by the 
upgraded expression of the two frequencies of the radial and 
vertical oscillations 



K = j-^ (V2 + Q + 2a Ji + CJ2) 



(66) 



4.2 Bifurcations from the disk and from the thin 
tube orbit 

We now start to illustrate the main body of results concern- 
ing the orbit structure as determined by the main bifurca- 
tions. 

Lowering the angular momentum (namely increasing 
the fictitious energy E) both normal modes may lose their 
stability through a 1:1 resonance. We denote the bifurcat- 
ing family as the inclined orbit in view of its natural in- 
terpretation as t he in phase (-i/) = 0) 1:1 r esonance of the 
two o scillations l|de Zeeuw fc Merrittlll983l : iBelmonte et al.l 
l2007h . The anti-phase (1/) — tt) 1:1 resonant loops never 
appear in these systems, at least as a stable family (see 
subsection I4.3|l . The inclined periodic orbits parents two 
families of inclined boxes that may arrive quite far from 
the equatorial plane both above and below the disk: th is 
phenomenon is called levitation (|Sridhar fc Tournal Il996l ). 
We recall that our inclin ed orbits have also been ref erred 
to as reflecte d banana bv lLees fc Schwarzschild (Il992l) an d 
lEvansI (Il994h and simply as banana by Hunter et al. ( 1998h : 
we p refer to leave this term as the stand ard denomina- 
tion jMiralda-Escude fc Schwarzschild|[l989l ) for the 2:1 res- 
onance. 

To describe the bifurcation of the inclined or- 
bit, the normal form is computed by a small 'detun- 
ing' (IContopoulos fc Moutsoulasiri966l : Ide Zeeuw fc MerrittI 
Il983h of the 1:1 resonance so that, from (|3ip with mi — 
m2 = 1, 



= qV2 + a = 1 + 5. 



(68) 



The normal form truncated to the first non zero resonant 
term is 



qK = Ji+J2 + SJi + q{ajl + 6 j| + cJi J2) (69) 



-|-dJi,72Cos[2(6li -02) 



K 



d= ^(l + «)(2-a). 



(70) 



(71) 



The rescaling by q lets to embody the detuned resonance 
(|68l) in the neatest form. Introducing the resonant combina- 
tions 

= 2(01 - 6*2) (72) 
and 

7^=Jl-J2, (73) 
the new Hamiltonian can be expressed in the reduced form 

k = £ + ^5{£+n)+A{£'^ + V?) + BE'Tl+ (74) 



and 



i(iS^ — TZ^){qc + dcos ip), 



(75) 



- + 2fe J2 + cJi 



(67) 



and by the orbit approximations (|51H52|I computed by 
means of the generating function (I60|l . These results are 
the basis f or attempting an accurate t racking of nor- 
mal modes (IContopoulos fc Seimenislll99d ) and box orbits 
l|Kent fc de Zeeuwlll99lC 



with 

A = |(a + 6), 
B = |(a-6) 
and 

£ = Ji + J 2 



(76) 
(77) 

(78) 
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Table 1. Critical value of the angular momentum for the bifur- 
cation of the inclined orbits from the disk orbit: the values in the 
third column are computed with Ea.l l86l l for general a and with 
ea. ll87l l for a = 0. The fourth and fift h columns list num erical 
predic t ions a nd t heir sources: jLces & S chwarzschiTdl lll992l) (LS); 
lEvanj il994) (E): iHunter et al.i (aQQa l (H); this work (T). 



g 


a 


Analytical 


Numerical 


Source 


0.4 


0.5 


0.27 


0.33 


H 


0.5 





0.30 


0.28 


T 


0.75 





0.42 


0.419 


LS 


0.8 





0.28 


0.26 


T 


0.8 


0.1 


0.24 


0.22 


H 


0.85 


-0.18 


0.23 


0.29 


E 



is the second integral as in (|45p . Considering the dynamics 
at a fixed value of we have that K defines a one-degree 
of freedom {tpjTZ) system with the following equations of 
motion 

i> = K'R^^5 + B£ + ^{A-{qc + dcosi;))TZ, (79) 



1Z = -K^ = ^ (£:^ - n') sin V. 



(80) 



The fixed points of this system give the periodic orbits of the 
original system. The pair of fixed points with TZ — ±£ corre- 
spond to the normal modes. The periodic orbits in "general 
position" (inclined and loops) are respectively given by the 
conditions ^ = and -0 — ±7r (which are the solution of 
TZ = 0) and by the corresponding solutions of ip — 0: 

S + 2BS 



TZ 

qc- 4A±d 
In view of (1731) and 



78l). the conditions 



< Ji, J2 sS £, 

applied to the solutions 
existence 



(81) 



(82) 



qc 



2b ±d 



\8T\l translate into the conditions of 



(83) 



and 

f > - 



2a 



The plus sign in front of d corresponds to the loops and gives 
a case that will be examined in the following subsection. The 
conditions with the minus sign in front of d correspond to 
the bifurca tion of the inclined. P roceeding along the lines 
followed bv lBelmonte et al] (|2007l ). it can be proven that, at 
the bifurcation, one of the normal modes suffers a stability- 
instability transition. The more common situation (for mod- 
els ranging from sensibly oblate to prolate) is that in which 
the 2;-axis becomes unstable and the inclined appears as a 
pitchfork bifurcation from the disk orbit. The passage to 
instability of the y-axis is possible only for strongly oblate 
models and gives rise to a pitchfork bifurcation from the 
thin tube. 

The bifurcation equations (I83I84|I determine critical val- 
ues of £ in terms of the parameters q,a. To make a quantita- 
tive prediction, we wont an expression for the corresponding 
critical angular momentum. The approach we have followed 



so far is altogether a perturbation approach truncated to 
the first non-trivial order. Therefore it is natural to look 
for expansions truncated to the first order in the detuning 
parameter. Taking into account the rescaling in K and the 
expressions (I631465|) and ((7T|, the first order expansions of 
the critical values of the fictitious energy (|19p are: 



f 12(2+a) , 

p ^ I 5(-2-a+Q^)' 

] ^^5, S>0. 



5<0, 



(85) 



The first solution corresponds to the bifurcation from the 
thin tube, the second one corresponds to the bifurcation 
from the disk. These are examples of series of the form (|55|l 
truncated to the first order. 

Using the relation between E and L established by (|19|l , 
we get the following expressions for the critical values of the 
angular momentum below which inclined orbits exist: 

2 + c. 



. _ 1 

i^crit — — p X ■ 



, 24a(q,/2 + a-l) 
^ 5{2-|-Q-c«^) 

, , 12Q(qy2+5'-l) 



q < 



q > 



/2+a ■ 



(86) 



It is also useful to write the limiting case of the logarithmic 
potential {a — 0): 



e TT!' 



75' 



q > 



(87) 



V2- 



A comparison with the outcome of numerical determi- 
nations of the bifurcation threshold allows us to evaluate the 
accuracy of these analytical predictions. In Table [1] the crit- 
ical value of the angular momentum for the bifurcation of 
the inclined orbits, computed with Ea. l|86p for general a and 
with ea. (l87|l for q = 0, are compare d with numerical data 
obtained either from published works (|Lees fc Schwarzschildl 
1 19921 : lEvan^l 19941 : iHunter et al.lll998h or by numerical com- 
putations made for the present paper. In this case, the bi- 
furcation has been detected tracing the instability thresh- 
old of the normal mode by means of the Floquet method 
iBender fc Orszadll978l ). 

The accuracy is particularly good when the model is 
close to the exact resonance. Overall, the discrepancy lin- 
early grows with detuning, as can be expected in this first 
order approach. The first two lines represent two strongly 
oblate models with the thin tube to become unstable: in the 
rather extreme case with q — 0.4, a = 0.5 the detuning is 
5 = —0.37 and the relative error in the prediction is 18%. 
In all other cases, the disk becomes unstable and the qual- 
ity of the prediction can be represented by the case with 
q — 0.8, a = 0.1 when the detuning is 5 = 0.16 and the rela- 
tive error in the prediction is 8%. We may guess a prediction 



-(g^/2T^-l) 



-^'critjtruc 2 

which can be used to further improve the accuracy of (|86|l. 

4.3 Return to stability of the thin tube 

Pitchfork or period-doubling bifurcation are usually followed 
by a second stability change when the second reson ant fam- 
ily appears l|Miralda-Escude fc Schwarzschildlll989l ). In this 
setting, this possibility occurs if the loops appear. Using the 
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Table 2. Critical value of the angular momentum for the bifur- 
cation of the inner tube orbits from the disk orbit: the values in 
the third column are computed with Eq. JTOSj. The sources are as 
in the caption to Table [T] 







L 


crit 






a 


Analytical 


Numerical 


Source 


1.1 


0.25 


0.51 


0.41 


E 


1.1 





0.50 


0.37 


T 



bifurcation equations (|83I84|I with the plus sign in front of 
d and the explicit expressions of (|63H65p . the inequality can 
be satisfied only with values of the parameters correspond- 
ing to rather extreme oblate models. This implies a negative 
value of the detuning and a critical fictitious energy 

Z + ct — 

We get the following expression for the critical value of the 
angular momentum below which loops bifurcate from the 
thin tube: 



1 



8Q(g\/2^ 



11 



2 + a- 



(90) 



This is again a pitchfork bifurcation so that the thin tube 
regains its stability. The loops are unstable and, lowering 
the angular momentum below the critical value, remain un- 
stable for every reasonable combinations of the parameters. 
At the same time, the inclined tend to occupy an even larger 
fraction of phase space. Referring again to the models in ta- 
ble [T] in the case with q — 0. 4, a — 0.5 th e critical value for 
the return to stability is 0.32 (|Evanslll994l ) and (|90p predicts 
0.18. 



4.4 Bifurcation of the inner thin tube 

A different phenomenon is that of the bifurcation of the in- 
ner thin tube. This happens only in prolate models and is 
due to a period doubling bifurcation from the equatorial orbit 
which becomes unstable at the bifurcation. In our terminol- 
ogy this is a banana orbit actually related to the 2:1 reso- 
nance. In the full three-dimensional problem, it gives rise to 
a thin tube that is always closer to the symmetry axis that 
the 'ou ter' thin tube ex amined above which remains always 
stable (|de Zeeuwlll985l ). 

In this case, the detuning, from (|30|) and (1311) with mi = 
2, 7712 — 1, is 



S = qV2 + Q - 2. 



(91) 



The normal form truncated to the first non zero resonant 
term is 

(92) 
(93) 
(94) 



K = 2Ji + J2 + SJi + DVJ1J2 cos(6li - 262) 
with 

a - 2 



D = 



^/2(2 4-0)1/4 2^2(2 -Fa)i/4' 
Introducing the resonant combinations 

ip = ei- 202 

and 



7^=Jl-2J2, (95) 
the new Hamiltonian can be expressed in the reduced form 

K = £ + l5{2£ + 7^) + -^V2£ + TZ(£ - 27^) cos ^, (96) 
5 a/S 5 

where now 

£ = 2Ji+J2. (97) 

The fixed points corresponding to the periodic orbits in 
general position are given by the solutions of 



D{7£ + 67^) - 2V5V2£TnS = 0, ^ = 0, (98) 
D{7£ + 67^) + 2V5V2£TTZS = 0, ip = -K. (99) 

Combining these with the existence conditions 

0s;Ji^£:/2, o^j2^£, (100) 

gives the critical value 

S = (101) 

and, taking into account the rescaling implicit in the normal 
form. 



(102) 



~_ 2{2 + a) 2 
^-(23^^- 

We now have an example of a series of the form (I55p in 
which the first coefficient vanishes and is therefore truncated 
to the second order. Using (|19|l the prediction of the value 
of the angular momentum below which the inner thin tube 
bifurcates from the disk is 

The pure logarithmic a = limit is given by 
LcHt=e-'*+*^^-^^'. (104) 

In Table [2] the critical value of the angular momentum for 
the bifurcation of the inner thin tube orbits, computed with 
Eq. p03p . are compared with numerical data. We may guess 
a prediction error 



: 2r = 2{qV2Ta-2y 



(105) 



^crit ,true 

which can be used to further improve the accuracy of (|103p . 

4.5 Bifurcations of higher-order resonant orbits 

As an application of the approach adopted bv lSanders et al.] 
(|2007l ) we compute the appearance of the pretzel as a 
very high-order (4:3) resonanc e. In the literarure it i s also 
denoted as a 'refiected' fish ijLees fc SchwarzschiTdl Il992l : 
Evans 1994). In spite of the high order, this family hap- 
pens to pay a relevant role in shaping the phase space of 
these systems, occupying a substantial fraction of surfaces 
of s ection for a wide range of parameters in oblate mod- 
els (iLees fc Schwarzschildi il992l : lEvand Il994l : iHunter et al] 
199% 

The implementation of the general procedure followed 
in all cases treated above would require, in the present in- 
stance, a normal form truncated to A'min = 7 in order to in- 
clude at least the first resonant terms. The algebraic manip- 
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ulators available nowadays have no problem in accomplish- 
ing this endeavour without excessive CPU times. However, 
this implies a huge number of terms which hinders a com- 
pletely general algebraic approach. The main troubles come 
from the solution of the equations for the fixed points of the 
reduced system, giving in turn the periodic orbits in general 
position. Being polynomials of degree one less than the nor- 
mal form in the action variables, these give rise to algebraic 
equations of which it is very difficult to write the solutions 
in a manageable way. A straighforward approach would be 
that of examining specific cases and numerically solving for 
their roots. However, this spoils this general method of all its 
appeal. The alternative is to give up for a detailed knowledge 
of the periodic orbits in general position and try an approx- 
imate location of their 'resonance manifold' together w ith 
its first order condition of existence jSanders et al.|[2007l l. 

We recall that, for a general reduced Hamiltonian of 
the form K{TZ,tp;£,q,a), the fixed points with ^ = and 
= ±7r correspond to two different orbit sets with two given 
values of TZ{0) and TZ{tv). They can be separately identi- 
fied only if the normal form contains at least the first res- 
onant term. The method adopted by iSanders et al. (2007), 
by truncating at a lower order term, only allows us to find 
an 'average' value of the two TZ coordinates together with 
an existence condition, a linear or quadratic equation in the 
simplest instances. However, it involves the detuning and 
the second approximate integral and therefore provides an 
estimate of the critical energy for the bifurcation. 

In this case the detuning, from (|30|l and pip with mi = 
4, 7712 ~ 3, is 



qV2r 



The truncated normal form is 



if = 4Ji + 3J2 + 35 Ji + 3q{aJt + bJi + CJ1J2) 



(106) 



(107) 



where the second order term is essentially the same as in 
the non-resonant normal form (|56p and the dots stand as a 
remainder of the possible necessity of continuing the expan- 
sion in case the quality of the approximation is too low. In 
this case one can simply chose 



TZ = 
and 



Ji 



(108) 



£^4:Ji+3J2. (109) 
The non constant terms of the reduced Hamiltonian are 



35TZ 



16b . \ ^2 /86 \ 
a + — —icj TZ^ + I — - c \ £TZ 



(110) 



The resonance manifold is determined by the condition tp 
K-jz = so that we have 



TZ 



275 - 8b£q + 3c£q 



2(9a+ 16fe- 12c)g' 
Combining this with the existence condition 
1 



7£ ^ TZf 



^ 0, 



(111) 



(112) 



the rescaling E — 3q£ and the expressions (|63H65p . we get 
the following second order approximation for the critical 'en- 
ergy' 



Table 3. Critical value of the angular momentum for the bifur- 
cation of the 4:3 pretzel (reflected fish). The sources are as in the 
caption to Table [l] 







L 


crit 




Q 


a 


Analytical 


Numerical 


Source 


0.75 





0.17 


0.173 


LS 


0.8 


-0.3 


0.10 


0.15 


H 


0.85 


-0.18 


0.18 


0.16 


E 


0.9 





0.40 


0.38 


T 



Ebm — — 



15 



v/2(2 + . 



-.6- 



—c 

2 + a 



(113) 



Using (|19|) . the prediction of the value of the angular mo- 
mentum below which the inner thin tube bifurcates from the 
disk is 



- f3ERM{q, a] 



(114) 



In Table [3] the critical value of the angular momentum 
for the bifurcation of the inclined orbits, computed with 
Ea. (|114p . are compared with numerical data. The accuracy 
reached is unexpectedly good in view of the rough nature 
of the approximating technique. This method could also be 
employed to investig ate the properties of resona nt boxlets 
of triaxial ellipsoids (|Zhaolll999l : lzhao et al.]|l999l '). 



5 DISCUSSION AND CONCLUSIONS 

Motion in the meridional plane of an axisymmetric system is 
a stumbling block and a case study in dynamical astronomy. 
Its investigation has a long history and has inspired the de- 
velo pment of many ideas and methods |Contopoulosl 19631: 
Henon fc Heiles|[l963 : lHorilll966l : lGustavsonlll966l : IVerhulstl 
Il979i ) . The development of techniques to investigate the reg- 
ular and chaotic structure of phase-space and of approxi- 
mating integrals of motion (a 'third integral' in this case, 
in addition to E and L) pay much to the activity on this 
topic. In particular, a systematic investigation of the reg- 
ular dynamics by constructing Ham iltoni an 'normal forms ' 
started with the work of iHoril (|l966l ) and iGustavsonI (|l966l ) 
and, suitably traduced in algorithms of normalization, is at 
the basis of many subsequent works including the present 
one. 

Here we have shown how a simple first-order truncation 
of the normal form is able not only to convey qualitative 
information on the phase-space structure but also quanti- 
tative predictions with concrete meaning for applications in 
galactic dynamics. In specific situations, especially in non- 
scalefree models, painstaking extended numerical simula- 
tions are usually required to even get a flavor of what is 
going on. A simple albeit rough analytical recipe can there- 
fore be quite useful to get an average picture. With this 
approach, in place of numerical simulations quite heavy in 
view of the large parameter space, simple analytical formulas 
are presented for the computation of the main bifurcation 
thresholds of the problem. For each scalefree potential with 
slope given by a and ellipticity q, a critical value of the angu- 
lar momentum of the full tridimensional system determines 
the appearance of some new orbit family and, possibly, a 
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change in the stabiUty nature of the 'parent' orbit. The ap- 
proach can clearly be generalized to non-scalefree models at 
the price of a fourth parameter in addition to L, a and q. 

Overall, these results make us confident that, in order 
to improve the accuracy or to explore the features of minor 
families and/or higher-order resonances, it can be worth the 
effort of constructing and investigating complete resonant 
normal forms. However, it is important to remark that the 
asymptotic nature of the results implies a trade-off between 
accuracy and extent of the phase-space region that is mim- 
icked by the normal for m. This region is largest at some 
optimal truncation order l|Efthvmiopoulos et al.l [20041 ) that 
depends on the nature of the system and the order of the 
resonance. The accuracy in the predictions of e.g. the bifur- 
cation thresholds can be further improved with respect to 
that at the optimal order but only at the expense of recon- 
struction of dynamics. The choice of the whole strategy is 
then determined by the needs of the problem at hands. 
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